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Abstract 

A new deterministic mathematical model for the transmission dynamics of Middle 
East Respiratory Syndrome Coronavirus (MERS-CoV) is proposed and fully analyzed. 
The presented model exhibits a unique endemic equilibrium and there is no infection free 
equilibrium due to constant influx of latent immigrants. An invasion threshold parameter is 
derived and stability analysis of the full model and its two special cases is carried out. The 
impact of quarantine and isolation is assessed via threshold analysis approach, while the 
impact of immigration on the disease prevalence is discussed. Indeed, we showed that MERS- 
CoV can be controlled by quick isolation or monitoring close contacts and quarantining of 
suspected latent immigrants. Further, numerical simulations of the model reveal that the 
disease can be contained if these preventive measures are combined with high reduction of 
immigration rate. 


Subject Classification: (2010) 34D23, 93C15, 93D20 
Keyivords: Middle East Respiratory Syndrome, Isolation, Quarantine. 


1. Introduction 

Middle East Respiratory Syndrome Coronavirus (MERS-CoV) belongs 
to the subfamily Coronavirinae in lineage C of the genus Betacoronavirus. 
MERS-CoV is one of the four Coronavirus groups of which Severe Acute 
Respiratory Syndrome (SARS-Cov) belongs [9]. In September 2012, The 
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World Health Organization (WHO), reported the First cases of MERS- 
CoV [26], a novel respiratory disease, initially localized to Middle East 
countries. It has a high potential for transmission to close contacts, as 
observed in communities with sporadic exposure, health care facilities, as 
well as families having contact with infected members [10]. Coronaviruses 
are wide spread in Bats around the world but can also be found in 
many other animal species as well [20]. The first index was reported in 
Kingdom of Saudi Arabia (KSA) on 13th June, 2012, when unknown 
coronavirus isolated from a 60yr old man admitted in hospital with a 7 - 
day history of acute pneumonia, fever, cough, expectoration, short breath 
and subsequent renal failure which later resulted in fatal outcome [28]. 
Shortly thereafter, on 14 th September, 2012, the United Kingdom Health 
Protection Agency (HPA) Imported Fever Service was notified of a case 
of unexplained virus, later named human coronavirus England 1. It has 
been isolated from another patient with history of severe respiratory 
illness in an intensive care unit (ICU) of London hospital. The index has 
been transferred from Qatar and had a history of travel to Saudi Arabia 
[3]. The origin of the disease was traced back to an earlier time - April, 
2012, when health care workers in an ICU of a hospital in Zarqa, Jordan 
were confirmed to have been contracted with the novel coronavirus [9]. 
From June, 2012 to 25 th September, 2013, WHO reported total of 133 cases 
with 55 deaths from 9 countries [20]. Up to 14 th March, 2014, there were 
477 globally reported cases and alerts including additional 7 countries 
[15]. Globally, from September 2012 to 17 th April, 2013, WHO has been 
informed of a total of 243 laboratory-confirmed cases of infection with 
MERS-CoV, including 93 deaths [25]. It has been observed that all cases 
outside Middle East are traced back to the region either by travel to or 
from the affected areas or by direct or indirect contact with others who 
have a travel history to Middle East [1]. 

The epidemiological data available suggest that the infection is 
primarily zoonotic in nature, with possible cross-species transmission 
betweenhumans and animals (Bats, Camels). Some of the reasons advanced 
are: Close phylogenetic similarity between viruses isolated from humans, 
bats and camels; presence of viral neutralizing antibodies in dromedary 
camels for the past 2 decades from (EGP, JO, KSA, OM, Canary Island), 
and none in other livestocks [1, 13, 16, 21, 22]; close association between 
camels and humans as source of food, sporting and as pets. However, 
most of the reported cases have no close contact with camels or bats . 
Therefore, the direction of transmission from humans to animals or vice 
varsa, host, reservoir e.t.c. remain unknown [1,16,17]. 
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The established clinical symptoms of MERS-CoV include fever, 
cough, shortness of breath, acute pneumonia, expectoration, hemoptysis 
often followed by renal failure [2, 20, 28]. Gastrointestinal symptoms like 
vomiting and diarrhea are also common [2]. The incubation period was 
estimated to be 5.2-5.5 days and the serial interval (the time between the 
successive onset of symptoms in a chain of transmission) was 7.6 days 
[2]. According to Cauchemez et al, [7], the estimated basic reproduction 
number obtained from all confirmed and probable of human cluster cases 
from Middle East (up to Aug 8, 2013) and the first 133 cases [20], were in 
the range 0.8-1.3 and 0.6 respectively. The implication of the former result 
is that self-sustaining transmission in the absence of control measures 
may occur. The average age of the first 133 reported MERS-CoV cases was 
52 years with male having higher case-fatality risk of 52% compared to 
women 24% [20]. 

Clinically, whenever there is an outbreak or re-emergence of new or 
an old infection in a community/country, in the absence of vaccine or 
treatment, isolation and quarantine of individuals with clinical symptoms 
or suspected to be exposed to the pathogens have proved to be effective 
ways of containing the disease [11, 13]. The World Health Organization 
is concerned, among others, with the gaps in understanding the degree 
of transmissibility between people, possibility of "super spreaders" and 
potential for sustainable person-to-person transmission and spread [10]. 
However, several mathematical epidemic models incorporating one of 
such control measures or both exist in the literature, see for instance, [11, 
18,23,24,27], For the 2015 epidemic outbreak of MERS-Cov in the Republic 
of Korea, two deterministic models for the transmission dynamics of 
MERS-Cov with and without control measures were developed in [27], 
The basic reproduction number, 1Z 0 was estimated to reach up to 4.422. 
The numerical analysis reveals that lack of self-protection sense and 
targeted control measures were the reasons of the outbreak spread quickly. 
However, it was reported that strengthening self-protection ability of 
susceptible and quickly isolating or monitoring close contacts are effective 
measures to control the disease. Furthermore, partial correlation analysis 
shows that the infectivity and proportion of the asymptomatic infected 
cases have much influence on the disease spread. In the same vain, another 
deterministic model was designed to analyze the MERS-CoV outbreak in 
the Republic of Korea [18]. The proposed model explicitly incorporates 
superspreading events and time-dependent transmission and isolation 
rates. The superspreaders are those who transmit the virus to more than 
20 patients and have underlying respiratory diseases with a severe cough. 
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The impact of the timing of control measures associated with a reduction 
of the transmission rate and diagnostic delays on the outbreak size and 
duration was assessed. It was observed that early interventions within 1 
week after the epidemic onset is a promising means to reduce the size and 
duration of the MERS-CoV epidemic. 

In a similar note, an SEIQJR - SARS model with both quarantine 
and isolation strategies was presented in [11]. The analysis of this model 
reveals that a reduction in the contact rate between susceptible and infected 
individuals by isolating the latter is a critically important strategy that 
can control SARS outbreaks with or without quarantine. It was reported 
that an optimal isolation is more effective than sub-optimal isolation 
and quarantine together. And that an optimal isolation combined with a 
highly effective screening programme at the points of entry would lead to 
a community-wide eradication of SARS. 

It is well known that migration and immigration greatly increase 
the spread of many infectious diseases at a regional, national and global 
scale [4, 12]. Simple models for disease transmission incorporating 
immigration of infective individuals was presented in [4]. The models 
exhibit a single endemic equilibrium that is asymptotically stable. It 
was revealed that there is no disease free equilibrium in the presence of 
immigration of infectives. For HIV transmission model a considerable 
reduction of infectives was suggested by screening and quarantining of 
infectives. While both immigration/migration terms were incorporated 
into all sub-population compartments, susceptible and infected, of two 
types of well-known heterogeneous epidemic models: multistage models 
and multi-group models for HIV/AIDS and other sexually transmitted 
diseases (STDs). It was shown that, the disease always becomes endemic 
in the population and tends to a unique asymptotically stable endemic 
equilibrium when migration or immigration into infected sub-population 
is present [12]. 

In this paper, a disease transmission model for MERS-CoV epidemic 
incorporating immigration of susceptible and latently infected individuals 
is proposed to assess the impacts of quarantine and isolation strategies in 
controlling the disease. As in the previous studies, the main result is that 
our model exhibit a unique asymptotically stable endemic equilibrium 
and there is no disease free equilibrium when immigration into latently 
infected sub-population is present. While such equilibria exist for the case 
when either all the immigrants are susceptible or the immigration rate is 
negligible. The impact of quarantine and isolation is assessed in relation 
to the relative infectiousness of quarantined and isolated individuals. 
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Further, eradication of MERS-CoV is feasible by reduction of immigration 
rate of latently infected individuals. The proposed model extends in some 
sense the aforementioned models in the literature as follows: 

(i) In the presented model, demographic features are incorporated 
and vital mathematical analysis is carried out, while these were not 
provided in MERS-CoV model proposed in [27]. 

(ii) In [18] there were only two sources of infection, here we consider 
four sources of infection. 

(iii) Although, immigration into both susceptible and infective 
compartments is common to [4, 12] and our model, in the latter 
we used saturated incidence formulation when the former used a 
simple mass action formulation. 

(iv) For the SARS model in [11], quarantined individuals are 
asymptomatic infectives who develop symptoms and then move to 
the isolated class. In our model, some quarantined individuals do 
not show symptoms and become susceptibles while others recovered 
by either treatment or acquiring natural immunity and move to 
recovered compartment. 

The remaining part of the paper is organized as follows. In section 
2, we describe the mathematical model. In section 3, rigorous analysis 
of the model is provided including its basic properties, existence and 
stability of the model equilibria, and the impact of quarantine, isolation 
and immigration rate are discussed. Numerical simulations supporting 
the theoretical results are provided. In section 4, we provide a concluding 
remarks. 

2. Model formulation 

The susceptible human population is generated via recruitment by 
birth (at a constant rate n) and immigration (at a rate aq); where 0 <q< 
1 account for the fraction of the inflow of migrants into the community 
(by air or road) who are susceptible, and by recovery of quarantined 
individuals at a rate <r(l - p); where 0 < p < 1 is a fraction of quarantined 
individuals tested with no clinical symptoms of the disease. This 
population is decreased following infection with MERS-CoV, which can 
be acquired via effective contact with infectious human (at a rate X), where 
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x _ Piu^+i+^Q+uJ) 

N 

Here the parameter /3 is the effective contact rate (contact capable 
of leading MERS-CoV infection). Furthermore, 0 < ri 1 < 1, 0 < r] 2 < 1 and 
0 < rj 3 < 1 are modification parameters accounting for the assumed 
reduction in infectiousness of individuals in the latent (L), quarantine (Q) 
and isolated (/) classes, in comparison to infectious individuals in I class. 
Natural death is assumed to occur in all human compartments at a rate q. 
Thus, the rate of change of the susceptible population is given by 

dS 

— = Tl + aq + a(l-p)Q-(A +ju)S. 

dt 

The population of individuals in the latent class is increased by 
immigration at the rate a( 1 - q) and by infection (at the rate A), and is 
decreased by the development of clinical symptoms (at a rate tJ, 
quarantine (at a rate r 2 ) and natural death, so that 

a(l-q) + AS-(r 1 +T 2 + / u)L. 


The population of infectious individuals with clinical symptoms 
of MERS-CoV in I class increases following the development of clinical 
symptoms by individuals in latent class (at the rate Tj). This population 
is decreased by isolation (at a rate y), recovery (at a rate SJ, MERS-CoV 
induced mortality (at a rate p) and natural death, this gives 

-77 = ^L-0 ' + v 1 + 8 l +/u)I. 


The population of quarantine individuals increases following the 
quarantine of individuals in the latent class (at the rate r 2 ). This population 
is decreased by recovery (at a rate a) and natural death, so that 


dQ 

dt 


r 2 L-(a + ju)Q. 


The population of individuals that are isolated (/(f)) is generated by 
the isolations of infectious individuals with clinical symptoms of MERS- 
CoV (at the rate y). It is decreased by recovery (at a rate <5 2 ), disease induced 
death (at a rate v 2 ) and natural death. Hence 
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Figure 1 

The flow diagram of model (1). The parameter, X = L + I +j2Q+'?3/) 


77 = ri -(v 2 +<5 2 + m)J- 

dt 


Finally, the recovered population is generated by the recovery of 
individuals in I, ] and Q classes at the rates <5,, S 2 and op, respectively, and 
decrease due to natural death. This gives 

^ = SJ + SJ + opQ-^R. 


In summary, the MERS-CoV transmission model is given by the 
following system of non-linear differential equations (a flow diagram of 
the model is given in Figure 1 and the associated variables and parameters 
are described in Table 1, respectively). 


dS 

dt 

dL 

dt 

<H_ 

dt 

dQ 

dt 


Fl + ac] + cr(l — p)Q — 


fiOpL + I + ^Q + iiJ) 
N 


S- juS, 


n \ /?('7i L + f + '7 2 Q + '7 3 /)c , u 

a(l-q) + - — -S-(p +t 2 +ju)L, 

t 2 L-(/ + v 1 +S 1 + h)1, 
t 2 L-(o + ju)Q, 


(1) 
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+J 2 +ju)J, 
at 


dR 


— = SI + SJ + crpQ-juR, 
dt 

with nonnegative initial conditions S(0) > 0,L(0) > 0,1(0) > 0,Q(0) > 0, 
J(0) > 0,R(0) > 0. 


Since the last compartment (R) does not appear in the first five 
equations of model (1), it is convenient to discard the last equation of 
system (1) for theoretical analysis. 

The equation for the total population which is obtained by adding all 
the equations of (1) is given by 


^ = n + a-//N(0-v 1 J(t)-v 2 7(t). 


( 2 ) 


3. Analysis of model (1) 

3.1 Basic properties 

The basic properties of the model can now be investigated. 

Lemma 1: The solution (S(t),L(t),I(t),Q(t),J(t)),R(t)) of model (2) with positive 
initial conditions, exists for all t > 0 and is unique. Furthermore, S(t) > 0, 
L(t) > 0,I(t) > 0,Q(t) > 0 and J(t ) > Ofor all t>0. 

Using Lemma 1 on equation (2), we have the following result. 

Lemma 2 : The biologically-feasible region of model (2) is 


T = {( S(t),L(t),I(t),Q(t),}(t),RmsR 6 + :S(t) + L(t ) 



which is positively-invariant and attracting. 


Proof: Adding the equations of the model (1) gives 


= n + a- //N(t) - Vjl(f) - v 2 J(t), 


so that 
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dN(t) 

dt 


n + a- /jN(t). 


It follows that dN ^ t ' > <0 if N(t) > . Thus, a standard comparison 

theorem can be used to show that 

N(t) < N(0)e“'" f + Htt^.[i_ g -^] 

M 

In particular, N(t) < if 7V(0) < n ^ a . Thus T is positively 

invariant. Furthermore, if N(t) > n+ ", then either the solution enters T 

n 

in finite time or N(t) approaches n+a and the infected variables L, I, Q 

and / approaches zero. Hence, T is attracting. That is, all solutions in R* 
eventually enter T. 


Table 1 

Description of parameters of model (1) 


Parameter 

Interpretation 

n 

Recruitment rate of susceptible humans by birth 

(] 

Fraction of immigrant humans 

a 

Recruitment rate of humans by immigration 

G 

Removal rate from quarantined 

p 

Transmission coefficient 

Uv V 2, ri 3 

Modification parameters for reduction in infectiousness of 
latent, quarantined & isolated individuals, respectively 

H 

Natural death rate in all classes 

T l> T 2 

Progression rates to infectious/quarantined classes, 
respectively 

r 

Isolation rate 

v„ V 2 

Disease-induced death rates for infectious and isolated 
individuals, respectively 

s v 4 

Recovery rates for infectious and isolated individuals , 
respectively 

V 

Fraction of quarantined humans 
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3.2 Existence and stability of equilibria 

It should be noted that model (1) has no disease free equilibrium (DFE) 

due to the constant inflow of newly latent individuals. Let £ = (S, L, I, Q, ]) 
be the equilibrium solution of system (1). Then setting the right hand side 
of (1) to zero we obtain 


S = 

L = 

i = 
Q = 
J = 


K 1 K 3 (U+aq) + aar 2 {l-p){l-q) 

(M + X)[r 2 p + (r, +ju)(ct + p)] + t 2 ct{ju + ip)' 
K 3 [A(Jl + a) + iia{l-q )] 

(M + + ( r i + /^)( CT + /*)] + r 2 <r(/z + h)' 

r 1 X 3 [i(Il + a) + pa(l-q)] 

K 2 {{p + A)[r 2 p + (r, + p)(a + p)\ + r 2 a{p + Ap)}‘ 
t 2 [A(J1 + a) + pa{\-q)\ 

(M + X)[t 2 ju + (r, + p)(a + p)] + r 2 a{p + ip)' 
t 1 /K 3 [i(tl + a) + jua(l - q)] 


(3) 


K 2 K 4 {{p + A)[r 2 p + (rj + p)(a + p)\ + t 2 cj{ju + Ap)} 

where K 1 =T 1 + r 2 +p,K 2 =y + v 1 + S 1 +p,K 3 =a + p, K 4 = v 2 + § 2 + p, 
f(i + /i 1 L + /j 2 Q + i iJ) 


/l = 


N 


and 


N = S + L + I + Q + J. 


(4) 


(5) 


Substituting (3) and (5) into equation (4) yields after algebraic 
manipulations (6). 

X{K 3 K 3 (II + aq) + aar 2 (1 - p)(l - q)} 


{p + A)[z 2 p + (t-j + p)(a + p)} + t 2 (t(p + Ap) 


+ OKU-PJ3) 


+<h 


K n 


+ t 2 (A-ii 2 P) + 




(6) 


k 2 k 4 


= 0 


o = 


[i(n + a) + pa{ 1 - q)] 


{p + A)[z 2 p + (z-j + p)(a + p)] + z 2 a(p + Ap) 


with 
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Thus, the endemic equilibria of system (1) correspond to the positive 
values of (3) obtained by using the positive solutions of (6). 

We define an invasion threshold by 

+ a)[r 1 K 3 (K i + r] 3 y) + K 2 K 4 (^K 3 + tj 2 t 2 )] 


%i K^KfU + aq) 

Thus, we obtain the following results. 


( 7 ) 


Theorem 1: Model ( 1 ) has: 

( i ) a unique endemic equilibrium for all values ofK i ifq e [0, 1) and a> 0, 

(ii) a unique endemic equilibrium ifTZ i > 1 and either q = 1 or a= 0, 

(iii) no endemic equilibrium ifTZ t < 1 and either q = 1 or a= 0. 

Proof : Re-arranging and simplifying equation (6) gives the following 
quadratic equation for X 

aX 2 +bX + c =0, (8) 


where 


a = (n + a)[R 2 R 4 (X 3 +r 2 ) + r 1 R 3 (K 4 +y)] / 
b = Y[K i (n + aq)(l-TZ l ) + a(l-q) J u{T 2 K 2 [K 4 a(l-p) + K lM )] 


+ K 3 K M T 1 +K i) + T 1 K 3 «>“]}/ 


3 4 ' ' 1 2 ' 13 

c = -af/./ 2 (1 - q)[K 2 K i (r/ 1 K 3 +r/ 2 r 2 ) + T 1 K 3 (K i +77 3 y)]. 


( 9 ) 


The coefficient a of equation (8) is always positive and c is always 
negative (noting that a > 0 and q e (0, 1)), whereas b is positive if 1Z i < 1. 
If TZ t > 1, b can also be positive or negative. Thus, model (1) has a unique 
endemic equilibrium for any sign of b by Descartes rule of sign. This 
proves (i). Next, if q = 1 or a = 0, then the coefficient ** and equation (8) 
reduces to 


aX + b =0, 
with either 

a = (n + a)[K^{K 3 +T 2 ) + T 1 K 3 {K t +y)l 
b = nfC(n + a)(l-7T), 

l =i 


( 10 ) 

( 11 ) 


or 
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a = n [K l K,(K 3 +T 2 ) + T 1 K 3 (K t + r)l (12) 

b = f\K n(i-n.), 

1 =1 

respectively. It is clear that, in each case the solution X = —~ of (10) is 
positive if > 1 and negative when lZ t < 1. Moreover, for TZ i = 1, the 
coefficient b of (10) is zero and so, the coefficient a of this equation must 
be zero since A cannot be zero. Hence, equation (10) has no solution 
whenever lZ j = 1. Consequently these prove (ii) and (iii). □ 


It follows from (7) that if either q = 1 or a = 0, then the invasion 
threshold, 1Z t reduces to 


n: = 

l 


Pl^K 3 (K 4 +ri 3 y) + K 2 K i (ii 1 K 3 +?i 2 T 2 )] 


^k 2 k 3 k 4 


(13) 


The uniqueness of the endemic equilibrium, £ allows us to 
numerically show its global stability as depicted in Figure 2, since the 
analytical proof is quite involved. 


3.3 Immigration cases 

As mentioned earlier, when all the immigrants are susceptible, 
q = 1 or the immigration rate is negligible (a = 0) the threshold parameter 
is as defined in (13). Furthermore, for each of such scenarios a disease free 
equilibrium exists. 

For the first case, in the absence of infection, model (1) with q = 1 
exhibits a locally asymptotically stable disease free equilibrium given by 


£ q = 
°o 


(S',L',r,Q',f) = 


n + a 
M 


,0,0,0,0 


In fact, the Jacobian matrix evaluated at £ 0 , denoted by 3„ is 


r 

~P 7 li 

p 

cr(l-p) 

-Pl 3 

0 

P 1 l\ - K 3 

-p 

Pi2 

Pi 3 

0 


-K 2 

0 

0 

0 

T 2 

0 

-*3 

0 


0 

r 

0 

-k 4 


(14) 







MERS-COV WITH LATENT IMMIGRANTS 


13 


o 



o o o o 

O O •T') 

<N 

uoii^[ndoj pappjuj 




o 



o o o o 

o o m 

<N ’— »— 1 

UOTJEjndOJ P9JD9JUJ 


IN 

d 


{{ 


o 


$ 

US' 

CO 

d 


'S' 


rH 

VA 

5-H 

o 

MH 





0> 

«+■* 

WD 

c 

'% 

o 

4= 


S 

z 


II 


0.71 and in (b) q = 0.52, a = 300, p = 3, 7?., = 9.55, H' = 8.49. 
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The eigenvalues of this matrix are: -fi, and the eigenvalues of the 
matrix 3, which is obtained by deleting the first/last rows and the first/ 
last columns of 3, • Then, the characteristic polynomial of this matrix is 

T 4 + AA 3 + BA 2 + CT + D = 0, (15) 


with 

A = fC 1+ tC 2 +fC 3 +fC 4 -/?; 7l 


B = 
C = 


k 3 k 4 + (K 3 + K 4 )(K 4 +K 3 -J3 ni )+K 2 (K 3 - Pnd- Ah + h7 2 ) 


K 1 K 2 K 4 


1 pT^Kt + nj) 


K 1 K 2 K 4 


+ K 1 K 3 (K 2 +K i ) 


1 A'/A+^c) 


*1*3 


+ K 2 K 3 K i -p(n 1 K 2 K i+ T 1 K 3 ) 

= K 3 K 2 K 4 ( i - 73 ,; )+K^K 3 (K 2 + tc 4 )(i - 73 ,;) + K 2 K 3 K 4 

-/3(? 1] K 2 K 4+ t ] K 3 ) 

D = K^K^il-n;). 


The Routh-Hurwitz criteria [19] for the local asymptotic stability in 
this case are A>0,C>0,D>0 and ABC-C 2 -A 2 D> 0. The first three 
conditions are satisfied if 73, < 1 which implies that 73, < 1 and 73 i2 < 1 
so that Pt 4 < K,f( 2 and /h;, < b,, respectively. The last condition holds 
since all terms remaining in the expansion after cancelation are positive. 
Thus the disease free equilibrium, £ 3 of model (1) with q = 1 is locally 
asymptotically stable when 73,* < 1. 

It is instructive to note that if 73* = 1, the coefficient D of the 
characteristic equation (15) is zero and so, J has a zero eigenvalue. This 
indicates the occurrence of a bifurcation when 73* = 1 which will be 

l 

investigated using Centre Manifold theory as described in Castillo-Chavez 
and Song [6, Theorem 4.1] in Subsection 3.4. 

The epidemiological interpretation of this result is that MERS-CoV 
will be eliminated if the initial sizes of the infected sub-populations of 
the model are in the basin of attraction of the DFE (£’). However, disease 
eradication is independent of the initial sizes if the DFE is globally 
asymptotically stable. The global stability result is established in the 
following theorem. 


Theorem 2 : The DFE (£„) of model (2) with q = 1 is globally-asymptotically 
stable (GAS) when 73,* < 2. 
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Proof: The proof is based on the approach in [1]. Using the notation in [5], 

for system (1) with q = 1,X = S,Z = (L,I,Q,J), and X*=j^i^,oJ. Then, 
the equation of uninfected is 


dX 

dt 


F(X, Z) = II + a + cr(l - p)Q -AS- /uS, 


and the infected subsystem is 


dZ 

dt 


G(X,Z) 


' AS-{t 1 + r 2 + ju)L 
T 1 L-(r + v 1 +5 1 + p)I 
r 2 L-(a + ju)Q 
jI-(v 2 +S 2 +ju)J , 


When L-I = Q = J = 0 (i.e., Z = 0), then F(X,0) = (II + a)-juS, and 


A = D G(X*,0) = 



P 


Pn 3 

h 

-k 2 

0 

0 

u 

0 


0 

1° 


0 

~k 4 


is obviously an M-matrix (a matrix whose off diagonal elements are 
nonnegative). 

Next, we have 

G(X, Z) = AZ - G(X, Z), which implies that 



s'). 

f s) _ 

f s') 

s V 

\ 

i 

1-— \+I 

{ N) 

1 - + r/,Q 

{ NJ 2 

i nJ + ">' 

I'-n) 



G(X,Z) = 


v 


0 

0 


It is clear that G(X, Z) > 0, since 0 < S < N. It is also clear that 

X' = n+ “,0j is a global asymptotic stable equilibrium of = F(X,0). 

Thus, by using a result in [5], we conclude that £„ is globally-asymptotically 
stable when 1Z* < 1. □ 


It follows from Theorem 1 that, model (1) with q = 1 has a unique 
endemic equilibrium if TZ* > 1. 
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Now, let and £ be the disease free and endemic equilibria of 
model (1), respectively with negligible immigration rate (a = 0), then we 
obtain the following results. 

Theorem 3 : The unique DFE, E“ of model (1) zvith a = 0 is globally- 
asymptotically stable iflZ* < 1, otherzuise unstable. 

Theorem 4 : The unique endemic equilibrium, £ a of model (1) zvith a = 0 is 
globally-asymptotically stable iflZ* > 1. 

Hence, the summary of the stability results of the model equilibria is 
given in Table 2. 


3.4 Bifurcation analysis 


Here we investigate the existence of bifurcation of model (1) using 
Centre Manifold theory since the Jacobian matrix J for the case when 
all the immigrants are susceptible has a simple zero eigenvalue when 
1Z* = 1 with other eigenvalues having negative real parts as presented in 
Subsection 3.3. In order to apply such a theory, we make the following 
change of variables S = x,, L = x 2 , 1 = x 3 , Q = x 4 , / = x 5 , so that we use the 
vector notation x = (x 1 ,x 2 ,x 3 ,x 4 ,x 5 ) T . Then system (1) with q = 1 can be 

written in the form = F(x) with / = (/i// 2 // 3 // 4 / f 5 Y> such that 


x 


2 


X 

X 


3 

4 


n + a + <t( 1 - p)x 4 


f(/) 1 x 2 +x 3 +y,x 4 +ri x 5 ) 

- — -—-- jux 

X , +x. +x, +x. +x c +x c 

1 Z 5 4 D o 


P(ll 1 X 2 +X 3 +ll 2 X 4 +ll 3 X b ) 
X 1 +X 2 +X 3 +X i +X b + X 6 

z 1 x 2 -(r+v 1 + S 1 +ju)x 3 , 
t 2 x 2 -(<t + ju)x 4 . 


X 1 -{t i +T 2 + fj)x 2 , 


yx 3 -{v 2 +S 2 +ju)x 5 . 


(16) 


The Jacobian matrix of system (16) at E 0 is the same with J as presented 
in (14). If ft is taken as a bifurcation parameter, we obtain from TZ* -1 that 

KKKK 

P = p =- 1 2 3 4 -. 

+ / I 3 r) + K 2 K 4 (/) 1 K 3 +t) 2 t 2 ) 

It can be verified that the Jacobian of (16) at ft = fS' has a right 
eigenvector associated with zero eigenvalue givenby u = {u l ,u 2 ,u 3 ,u 4 ,u 5 ) T , 
with 
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r 




G « s G ij yr 
-n + — + a(l-p) - 5 — 

1 V ' V V V 

K 2 ^2^4 J 


U 2 ,U 2 =u 2 > 0,U 3 = T 7 “M 2 , 
K 2 


G , y 

u. = —M 0/ and = —. 

4 x 3 2 5 k 4 


(17) 


Similarly, the left eigenvalue of J ? (f 0 ) associated with zero eigenvalue 
at p = P is given by v = (v 1 ,v 2 ,v 3 ,v 4 ,v 5 ) T , where 

v i =°' v 2 =v 2 >0,v 3 =-J^—(y, h -K 4 )v 2 ,v 4 =^v 2 ,v 5 =^v 2 . (18) 

R 2 R 4 k 3 k 4 

To establish the existence of bifurcation, we need to compute a and b 
as follows. Let F k be the k‘ h component of F, then as defined in [6] 


5 Q*p 

a = V v.u.u. - —(f) 

^ * ' ’dx .x 0 




i 


8 F, 


(19) 


b = V v u. -—(f ). 

ch.x . °' 


Now, for system (16) the associated non zero partial derivatives of F 
at the disease free equilibrium are as follows: 

8 % _ 2/?V _ /?>(77 1 + 1 ) 

( 20 ) 




IT + a 5x 2 x 3 


Il + a 


8 F 2 _ jui^+iiJ 8 F 2 0 Mfa+ils) 


8x 2 x 4 


n + a dx 2 x 5 


n + a 


It follows from (19) and (20) that 
Pfiu 2 v 2 


a = — 


Tl + a 


2ntl , h( ? 7i+ 1 ) , G(G + n 2 ) , r(G +/ 7 3 ) 


fC 




K. 


( 21 ) 


which is negative, and for the sign of b, it is associated with the following 
non zero partial derivatives of F 


8 2 F 8 2 F „ 8 2 F 8 2 F 

-— =-'7,, -— = -1, -— = - 77 ,, -— = 

8x 2 8p 1 5x 3 5/7 5x 4 5/7 2 8x b 8p 

8 2 F 0 2 F, „ d 2 F 0 2 F, 

5x 2 0yd 1 5x 3 3/7 8x 4 8p 2 dx 5 dp 


(22) 
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Table 2 

Summary of the stability results of the full model and its two special cases. 


Equilibrium 

Invasion threshold 

Stability 

£ 

TZ. f 1 

GAS 

% 

R*<1 

LAS/GAS 

% 

n'< i 

LAS/GAS 

i 

r*>i 

GAS 

£ a 

r*>i 

GAS 


It follows from the expression of b in (19) and the derivatives in (22) 

that 


r 


b = 


\ 


v K 2 K S 


v 2 > 0. 


(23) 


We obtain from the computations of a and b above that a is negative 
and b is positive. Thus, using Theorem 4.1 (iv) in [6], we establish the 
following result 


Theorem 5 : The system of equations (2) with q = 1 exhibits a transcritical 
bifurcation at £ 0 if the bifurcation parameter P' = or 

equivalently 1Z] = 1. 


3.5 The impact of quarantine/isolation/immigration 

In this section, the impact of quarantine and isolation on the 
transmission dynamics of MERS-CoV is investigated via a threshold 
analysis approach on the invasion threshold parameter, TZ i . This is carried 
out in terms of the parameters associated with quarantine of latent 
immigrants t 2 and the isolation of individuals with disease symptoms y. 
More precisely, the threshold TZ j is considered as a function of t 2 and y. 
For the case of the quarantine of latently infected individuals, the rate of 
change of 7?. : with respect to t 2 is as follows: 
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_ Pli P[K 2 K 4 (p 1 K 3 + , i 2 T 2 ) + T 1 K 3 (K i +i] 3 y)] 


dr 2 K,K 2 




Let us define from (24) that 

K^K.-r,) 
r]] K 2 K 3 K 4+ T,K 3 (K 4+ ri 3 r) 


n 2 = 


K 2 K 4 (t 1+M ) 
since K 1 =t 1 +t 2 + /u. Therefore, 


dll. 

L sg 0 if and only if // 2 ^ ? 7 2 


(24) 


(25) 


Thus, this threshold parameter will be a decreasing function of the 
quarantine parameter t 2 , resulting in a decrease of disease burden (new 
infections and mortality) if the relative infectiousness of quarantine 
individuals ri 2 does not exceed its associated threshold 77 *. By contrast, if 
r\ 2 < 77 *, then the use of quarantine of latent immigrants will increase the 
threshold parameter (Jl) and so, increase disease burden. In this case, the 
use of quarantine is detrimental to a community. The result is summarized 
as follows. 


Proposition 1: The use of quarantine of the latent immigrants will have positive 
(negative) impact in a community if r] 2 ;g 77 *. 

Similarly, the rate of change of H i with respect to isolation parameter 
yis given by 

dfc = +ii 2 t 2 ) + t 1 ii ? K 3 \ 

dy K,K 2 K 3 K 4 

p[K 2 K 4 (ii l K 3 +r l 2 r 2 )+T 1 K 3 (K i +T l 3 y)] 

v ir- is v ' ' 


We obtain from (26) that 

K 4 S 2 +V 2 +1 1 

Is = T-= TT-, 

b 1 +v l + /a b 1 +v 1 +11 


(27) 


so that 


- ^ 0 if and only if 
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It follows that, 7 Zj will be a decreasing function of the isolation 
parameter y and hence decrease disease burden if i] 3 < if. It will be an 
increasing function of y, on one hand, when >i 3 > if . In such a case, the use 
of isolation will have a detrimental impact in a community. This result is 
summarized in the following proposition. 

Proposition 2 : The use of isolation of infectious individuals zvill have positive 
(negative) impact in a community if p 3 

These results (Propositions 1 and 2) are illustrated numerically in 
section 3.5.1 below. 




(b) (c) 


Figure 3 

The simulations showing the impact of immigration on the number of infect¬ 
ed individuals using parameter values as in Table (3). (a) The infected humans 
increases with increasing influx of immigrants, (b) Indicates that the infected 
humans increases to a certain level, (c) Reveals that there is a critical immigra¬ 
tion rate above which the infected humans decreases. 
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It worthies mentioning here that, the critical parameters rj\ and rj\ as 
defined in (25) and (27), respectively, remain the same for the case when 
either all immigrants are susceptibles (q = 1) or the immigration rate is 
negligible, a = 0. Hence, the results in Propositions 1 and 2 also hold true 
for such two cases. 

To explore the impact of influx of immigrants into the susceptible and 
latent compartments on the dynamics of MERS-CoV, we use numerical 
experiment to vary 

different values of the immigration rate (a), while keeping all other 
model parameters fixed. In fact, the numerical simulations reveal that 
the number of infected individuals increases with increasing values of a. 
As depicted in Figure 3(a), when the influx increases from 0 to 700, the 
infected individuals increases from 30 to 90 as well. It also shows that 
the prevalence is very low (ranges between 10 to 14) and the steady state 
value I is also very small, whereas, for 0 < a < 100, it ranges between 32 to 
42 about, and if a = 700 then 1(200) is less than 100. Moreover, Figure 3(b) 
shows more clearly that the number of infected individuals increases with 
increasing value of a, while Figure 3(c) indicates that, there might be a 
critical value of a above which the disease prevalence decreases due to the 
positive impact of quarantine and isolation (Propositions 1 and 2) so that 
the disease will die out with time. However, eradication of MERS-CoV 
in a community is feasible if quarantine and isolation is combined with a 
high reduction of immigration rate. 




(a) Quarantine 


(b) Isolation 


Figure 4 

The simulations showing the impact of quarantine and isolation in controlling 
the infection with parameter values as in Table 3. (a) Quarantine with critical 
value, i] 2 = 0.3118. (b) Isolation with critical parameter 77 * = 0.1087 
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Table 3 

Parameter values for the model (1) 


Parameter 

Nominal value 

Reference 

n 

400 per day 

Assumed 

P 

0.25 per day 

Assumed 

a 

> 0 per day 

Assumed 

b 

0.0000374 per day 

[11] 

b, 

(0.001, 0.65) 

Assumed 

<i 2 

(0.013, 0.78) 

Assumed 

% 

(0.002, 0.65) 

Assumed 

Y 

(0.00138, 0.05) per day 

[2,8] 

O 

(0.00138, 0.025) per day 

[2, 8] 

h 

(0.0018, 0.0028)per day 

[2, 25] 

v 2 

(0.0012, 0.0018)per day 

[2, 8, 25] 

h 

(0.0018, 0.0028)per day 

[2] 

h 

(0.0012, 0.0018)per day 

[2, 8] 

<5, 

(0.0337, 0.26) per day 

[2] 

4 

(0.26, 0.5) per day 

[2] 

P 

[0,1] 

Assumed 

q 

[0,1] 

Assumed 


3.5.1 Numerical simulations 

The numerical simulations in Figure 4(a) is depicting the result of 
Preposition 1 in which the invasion threshold parameter lZ jr is a decreasing 
function of the quarantine parameter t 2 , whenever the modification 
parameter for quarantine is less than a critical value rf 2 =0.3118. Thus, 
for values of r) 2 = 0.017, 0.25, less than if 2 , the cumulative new cases of 
infection decreases from about 27,000 to 16,500 individuals. However, for 
rj 2 = 0.5, 0.75, greater than the critical value, the cumulative new cases 
increases from about 28,500 to 38,500 individuals. 

TZ t , is a decreasing function of the isolation parameter y whenever 
the modification parameter for quarantine is less than a critical value 
77*3 = 0.1087. Thus, for values of i] 2 = 0.001 to 0.017, less than T) v the 
cumulative new cases of infection decreases from about 18,500 to 16,500 
individuals. However, for rj 3 = 0.25 to 0.55, greater than the critical 
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value, the cumulative new cases increases from about 26,500 to 32,500 
individuals. 

4. Discussion and Conclusion 

A new deterministic model for the spread of the Middle East 
Respiratory Syndrome coronavirus in a community with quarantine and 
isolation as control strategies is formulated and rigourously analyzed. We 
notice that the presented model does not have disease free equilibrium 
point due to constant in flux of the latent immigrants and so, there is an 
endemic persistence of MERS-CoV depending on the invasion threshold 
parameter 7?.„ immigration rate as well as the epidemiological status 
of the immigrants. In fact, if the immigration rate a > 0 and fraction of 
susceptible immigrants q e [0,1), then the disease invades the population 
for all values of 1Z r On the other hand, endemic persistence is only 
possible if either immigration rate is negligible or when all immigrants 
are susceptible provided 71* > 1. Stability analysis is carried out for the full 
model and these two immigration scenarios. We showed numerically that 
this equilibrium point is a global attractor when TV V b see Figure 2. For 
the case when all the immigrants are susceptible, we proved the global 
asymptotic stability of the associated disease free equilibrium (£„) if the 
threshold parameter 71* < 1. Such a stability result hold true for the case 
when the immigration rate is negligible, a = 0 and the result follows easily 
from the first scenario. The existence and stability of endemic equilibrium 
for these scenarios follows from that of the full model. 

In order to complete our investigation, we assess the impact of 
quarantine and isolation via a threshold analysis approach. Indeed, we 
showed that quarantine of latent immigrants will decrease the burden of 
the epidemics (new infections), if the rate of the relative infectiousness of 
quarantined individuals does not exceed certain threshold quantity, if 2 . 
In a similar note, an isolation of infected individuals using intensive care 
units of hospitals, will have positive impact on the burden of the epidemics 
(new infections and mortality due to disease) by reducing the threshold 
parameter, if the relative infectiousness of isolated individuals does not 
exceed certain threshold quantity, rf 3 . These results are summarized in 
Propositions 1 and 2 and numerically supported by Figure 4. 

Furthermore, numerical simulations of the model reveal that the 
influx of susceptible and latent immigrants increase the number of 
infected individuals. In fact, one can observe from Figure 3(a) that if the 
influx of these immigrants increases from 0 to 700, the number of infected 
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individuals rises from 30 to 90. Despite such an effect of the immigration 
rate (a), the positive impact of isolation and quarantine drastically reduces 
the disease prevalence. As shown in Figure 3(a), the prevalence is very 
low (ranges between 10 to 14) and the steady state I is also very small. 
While for 0 < a < 100 it ranges between 32 to 42 about, and if a = 700 then 
1(200) is less than 100. Moreover, Figure 3(b) shows more clearly that the 
number of infected individuals increases as the immigration rate increases, 
whereas. Figure 3(c) indicates that, there might be a critical value of a 
above which the number of infectious humans decreases perhaps, due to 
the positive impact of quarantine and isolation (Propositions 1 and 2) so 
that the disease will die out with time. It is worthy of mention that, most 
of the parameter values used are obtained from some related literature 
as presented in Table 3. The ranges of such parameter values are chosen 
from minimum and maximum as presented in such existing works for the 
sake of consistency. While we assumed the values of the parameters ( p , 
q and q, for i = 1, 2,3) in the interval [0, 1] since each of such parameters 
is either a fraction of humans or a fraction that account for the assumed 
reduction in infectiousness of individuals in the latent, quarantined and 
isolated compartments. Moreover, the values of the recruitment rate II 
and the effective contact rate f3 are assumed in relation to their values used 
for SARS model in [11]. 

In conclusion, we showed using a deterministic model that MERS- 
CoV can be controlled by quick isolation or monitoring close contacts 
with the patients infected by MERS-CoV and quarantine of suspected 
latently infected individuals. This is similar to the finding in [27]. In 
addition, reducing the influx of immigrants plays a significant role in 
decreasing the number of infected individuals when the recruitment rate 
of immigrants is below a certain critical value. This is a value of a below 
which the prevalence is high. Flowever, community-wide eradication of 
MERS-CoV is feasible when quarantine and isolation are combined with 
high reduction of the immigration rate into a community. 
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